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O 
^J , Abstract 

We study a ferromagnetic Ising model on random graphs with a power-law degree distribution and 

compute the thermodynamic limit of the pressure when the mean degree is finite (degree exponent 

t > 2), for which the random graph has a tree-like structure. For this, we closely follow the analysis by 

Dembo and Montanari in [7J which assumes finite variance degrees (t > 3), adapting it when necessary 

Oh ' and also simplifying it when possible. Our results also apply in cases where the degree distribution 

r^ does not obey a power law. 

We further identify the thermodynamic limits of various physical quantities, such as the magneti- 
zation and the internal energy. 



(Nj '■ 1 Introduction and results 

>: 

In this article we study the behavior of the Ising model on complex networks. There are many real- 

iy-s [ world examples of complex networks. In [18], Newman divided such networks into four categories: social, 
information, technological and biological networks. There has been much interest in the functionality of 

l/~) ■ such networks in recent years [TJ \TE[ 122] . The Ising model is a paradigm model in statistical physics for 

O | cooperative behavior pjjl [20] . 

De Sanctis and Guerra studied this model on Erdos-Renyi random graphs in the high and zero tem- 
perature regime [21]. In [7J, Dembo and Montanari study a ferromagnetic Ising model on locally tree-like 
graphs, where they assume that the degree distribution of the graph has finite variance. The Ising model 

KS \ on the /c-regular graph where there is no external magnetic field is studied in more detail in [17]. In 
C$ ' this paper, the Gibbs measures are studied and it is proved that they converge to a symmetric linear 
combination of the plus and the minus Gibbs measure, while other Gibbs measures (of which there are 
uncountably many) are not seen. 

Many real-life networks are reported to have an infinite variance degree distribution (see e.g. |18j) 
and, therefore, it is interesting to generalize the analysis of the Ising model on random graphs to this 
setting. In this article we shall extend and simplify the analysis in [7] to the case where the variances 
of the degrees are infinite, but their means remain finite. In particular, we shall prove that the explicit 
expression for the pressure found in [7J remains valid in the case of infinite variance degrees. 

This research fits into a general effort to study the relation of the topology of networks and the 
behavior of processes on them. An overview of results by physicists can for example be found in |10j . 
Also mathematically rigorous results for processes on power-law random graphs were published recently, 
for example for the contact process [6] and first passage percolation [3]. 
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In this section we will first define the model and then state our main results. Furthermore we will 
discuss these results and give an overview of the proof. The remainder of the proof can be found in the 
subsequent sections. 

1.1 Model definition 

We start by defining Ising models on finite graphs. Consider a random graph sequence {G n }n>ii where 
G n = (V n ,E n ), with vertex set V n = [n] = {1, . . . ,n} and some random edge set E n . To each vertex 
i € [n] we assign an Ising spin <7j = ±1. A configuration of spins is denoted by a = {oi : i £ [n]}. The 
Ising model on G n is then defined by the Boltzmann distribution 

^ a "> = 7 (R m exp {P J2 ai<J 3 + /Z Bi(Ti ) ■ ( L1 ) 

ni ^'- j ( d,J)eE n ie[n] J 

Here, j3 > is the inverse temperature and B = {Bi : i 6 [n]} £ R n is the vector of external magnetic 
fields. We will write B instead of B_ for a uniform external field, i.e., Bi = B for all i G [n]. The partition 
function Z n ((3,B_) is the normalization factor: 

z n (p,B)= Yl ex p|^ Yl wj + Y, Bi<Ti \- ( L2 ) 

<T6{-1,+1}" [ (i,J)eE n i£[n] J 

We let ( • ) denote the expectation with respect to the Ising measure /j,, i.e., for every bounded function 
/:{-l,+l}»-+R, 

(f^), n = E f(°)»n(<7)- (1-3) 

ffS{-l,+l}" 

The main quantity we shall study is the pressure per particle, which is defined as 

MP, B) = - log Z n (/3,B), (1.4) 

n 

in the thermodynamic limit of n — > oo. 

We shall assume that the graph sequence {G n } n >i is locally like a homogeneous random tree, uniformly 
sparse and has a degree distribution with strongly finite mean. We make these assumptions precise below, 
but we shall first introduce some notation. 

For a probability distribution over the non-negative integers P = {Pk : k > 0} we define its size-biased 

law p = {pk ■ k > 0} by 

(fc + l)Pfe+i n ,, 

Pk = = , (1.5) 

where P = J2k>o ^Pk is the expected value of P. Similarly, we write ~p = J2k>o ^Pk f° r the expected 
value of p. The random rooted tree T(P, p, £) is a branching process with £ generations, where the root 
offspring has distribution P and the vertices in each next generation have offsprings that are independent 
and identically distributed (i.i.d.) with distribution p. We write P for the law of T(P, p, oo) and write 
T(p, £) when the offspring at the root also has distribution p. 

We write that an event A holds almost surely (a.s.) if ¥[A] = 1. The ball of radius r around vertex i, 
Bi(r), is defined as the graph induced by the vertices at graph distance at most r from vertex i. For two 
rooted trees 71 and 72, we write that 71 — 72, when there exists a bijective map from the vertices of 71 
to those of 72 that preserves the adjacency relations. 



Definition 1.1 (Local convergence to homogeneous trees). Let P n denote the law induced on the ball 
Bi(t) in G n centered at a uniformly chosen vertex i £ [n]. We say that the graph sequence {G n } n >i is 
locally tree-like with asymptotic degree distribution P when, for any rooted tree T with t generations, we 
have that, a.s., 

lim F n [Bi(t) ~ 71 = ¥[T(P,p,t) ~ 71. (1.6) 

n— too 

Note that this implies that the degree of a uniformly chosen vertex of the graph has asymptotic law 
P. In [7], it is assumed that the asymptotic degree distribution P has finite variance. This is not a 
necessary condition and we shall prove that it is sufficient to assume that the degree distribution has a 
finite (1 + e)-th moment for some e > 0: 

Definition 1.2 (Strongly finite mean degree distribution). We say that the degree distribution P has 
strongly finite mean when there exist constants r > 2 and c > such that 



J^Pi^ck-^-V. (i.7) 



i=k 



For technical reasons, we will assume, without loss of generality, that r £ (2, 3) in the rest of the 
paper. Note that all distributions P where 

oo 

J2 P i = ck ~ (T ~ 1)L (k), (1-8) 

i=k 

for c > 0, r > 2 and some slowly varying function L(k), have strongly finite mean, because by Potter's 
theorem ( |11| Lemma 2, p. 277]) any slowly varying function L(k) can be bounded above and below by an 
arbitrary small power of k. Also distributions which have a lighter tail than a power law, e.g. the Poisson 
distribution, have strongly finite mean. 

Definition 1.3 (Uniform sparsity). We say that the graph sequence {G n } n >i is uniformly sparse when, 



a.s., 



lim limsup- } Dili D>l \ = 0, (1.9) 



i^-oo n-S>oo n . 



where Di is the degree of vertex i in G n and \_\ denotes the indicator of the event A. 

An immediate consequence of the local convergence and the uniform sparsity condition is, that, a.s., 



I'-l - lim J_ y Ykt w ._ k} = - lim lim fvfc E * g[w]1{A=fc} + ly D1 

n^oo 2n ^ ^^ t^.-*l o ^->oon->oo ^ n n ^^ ' ^' 



lim .... .. , 

n^oo n n->oo 2n *— ' *— ' l 2 £->oo«->oo 1 ^^ n 77, 

iefnl fc=l \fc=l iefnl 



>n 



1 ^ _1 _ 

= - lim V fcP fc = P/2 < oo. (1.10) 

2 £— too *—* 

k=\ 

1.2 Main results 

We first investigate the thermodynamic limit of the pressure: 

Theorem 1.4 (Thermodynamic limit of the pressure). Assume that the random graph sequence {G n } n >\ 
is locally tree-like with asymptotic degree distribution P, where P has strongly finite mean, and is uniformly 



sparse. Then, for all < (3 < oo and all B € R, the thermodynamic limit of the pressure exists, a.s., and 
equals 

lim il> n (p,B) = <p(J3,B), (1.11) 

n— >oo 

w/iere, /or P < 0, p(/3, B) = M -B), <p(J3, 0) = lim B40 ¥>(A 5) and, for B > 0, 

P P 

</?(/?, P) =— logcosh(/3) - — E[log(l + tanh(/3)tanh(/ii)tanh(/i 2 ))] 



+ E 



log ( e B Y[{1 + tanh(/3) tanh(/i;)} + e~ B ]J{1 - tanh(/3) tanh(/i,;)} 



i=l i=l / 

u>/iere 

(%) L has distribution P; 
(ii) {hi}i>\ are i.i.d. copies of the fixed point h* = h*(/3,B) of the distributional recursion 

Kt 



(1.12) 



h (t+i) £ B + ^atanh(tanh(/3)tanh(/if ) )), (1.13) 

j=l 

where h^' = B, {Kt}t>\> are i-i-d. with distribution p and {h\ }»>i are i.i.d. copies of h^> inde- 
pendent of Kt; 

(Hi) L and {/ij}j>i are independent. 

The quantity ip(f3,B) can be seen as the infinite volume pressure of the Ising model on the random 
Bethe lattice, where every vertex has degree distributed as P (cf. |3j where the Ising model on the regular 
Bethe lattice is studied). 

Various thermodynamic quantities can be computed by taking the proper derivative of the function 
</?(/?, B) as we shall show in the next theorem. 

Theorem 1.5 (Thermodynamic quantities). Assume that the random graph sequence {G n } n >i is locally 
tree-like with asymptotic degree distribution P, where P has strongly finite mean, and is uniformly sparse. 
Then, for all (3 > and B ^ 0, each of the following statements holds a.s.: 

(a) Magnetization. Let M n (f3,B) = — EieM (°*)u ^ e the magnetization per vertex. Then, its ther- 

modynamic limit exists and is given by 

Af (0, B) = lim MM B) = — <p(j3, B). (1.14) 

n->oo OB 

(b) Internal energy. Let U n {j3,B) = —— E(j „•)<=£ ( a i a j) u be the internal energy per vertex. Then, its 

thermodynamic limit exists and is given by 

U(P, B) = lim UM, B) = ~M B). (1.15) 

n-*oo Op 

(c) Susceptibility. Let X n(P,B) = £Eije[n] (V^iU " (^n (^J = ^wi^ B ) be the suscepti- 

bility. Then, its thermodynamic limit exists and is given by 

X (P, B) = lim Xn (/3, B) = J^<p(p; B). (1.16) 



The limits above hold for j3 > and B ^ 0. Prom the physics literature (see e.g. [91 [16]) we expect 
that this model has a ferromagnetic phase transition at /3 C = atanh(l/p), i.e., the susceptibility becomes 
infinite at /3 = j3 c in B = 0. For (3 < j3 c the functions above are continuous in B for all B and thus the 
limits above also hold in this regime. 

Another physical quantity studied in the physics literature is the specific heat, 



C n {^B) = -0 



, dU n 

d(3 



(1.17) 



? 2 d 2 



Unfortunately, we were not able to prove that this converges to (3 -^gjip((3,B), because we do not have 
convexity or concavity of the internal energy in j3. We expect, however, that this limit also holds. 

Taking the derivatives of Theorem 11.51 we can also give explicit expressions for the magnetization and 
internal energy which have a physical interpretation: 

Corollary 1.6 (Explicit expressions for thermodynamic quantities). Assume that the graph sequence 
{G n } n >i is locally tree-like with asymptotic degree distribution P, where P has strongly finite mean, and 
is uniformly sparse. Then, for all /3 > and B £ M., each of the following statements holds a.s.: 

(a) Magnetization. Let vl+i be the random Ising measure on a tree with L + 1 vertices (one root and 
L leaves), where L has distribution P, defined by 



1 



"i+i" 



Z L+ i(P,h* 



■ exp 



L L \ 

j3 ^ a Q°~i + B(J o + ^2 hi(Ti f 



(1.18) 



where {/ij}j>i are i.i.d. copies of h* , independent of L. Then, the thermodynamic limit of the 
magnetization per vertex is given by 



vl+\ 



M(P,B)=E[(a ) 

where the expectation is taken over L and {/i«}j>i. More explicitly, 

/ l 



M(/3,B) =E 



tanh I B + > atanh(tanh(/3) tanh(/ij)) 



8=1 



(1.19) 



(1.20) 



(b) Internal energy. Let v' 2 be the random Ising measure on one edge, defined by 

1 



^) 



Z 2 (P,hi,h 2 ) 



exp {/3<7i<72 + hiai + /12C2} , 



(1.21) 



where h\ and h 2 are i. i. d. copies of h* . Then the thermodynamic limit of the internal energy per 
vertex is given by 



U(f3,B) = ~E 



\0-\u2) 



where the expectation is taken over h\ and h%. More explicitly, 

tanh(/3) + tanh(/ti) tanh(/i2) 



U{0,B) 



2 



1 + tanh(/3) tanh(/ii) tanh(/i2) 



(1.22) 



(1.23) 



Note that the magnetization and internal energy are local observables, i.e., they are spin or edge 
variables averaged out over the graph. This is not true for the susceptibility, which is an average over 
pairs of spins, and hence we were not able to give an explicit expression for this quantity. 



1.3 Discussion 

We study the Ising model on a random graph, which gives rise to a model with double randomness. Still, 
in the thermodynamic limit, the pressure is essentially deterministic. This is possible, because it suffices 
to study the Ising model on the local neighborhood of a uniformly chosen vertex. This local neighborhood 
converges by our assumptions to the tree T(P, p, oo), and it thus suffices to study the Ising model on this 
limiting object. An analysis of this kind is therefore known as the objective method introduced by Aldous 
and Steele in [2J. 

The assumption that the graph converges locally to a homogeneous tree holds in a wide range of 
random graph models, among which the configuration model and the Erdos-Renyi random graph, as we 
now explain. 

In the configuration model, a random graph G n is constructed as follows. Let {Di}f =1 be a sequence 
of i.i.d. random variables with a certain degree distribution P. Let vertex i £ [n] be a vertex with Di 
half-edges, also called stubs, attached to it, i.e., vertex i has degree D^. Let L n = Y27=i^i De ^ ne total 
degree, which we assume to be even in order to be able to construct a graph. When L n is odd we will 
increase the degree of D n by 1. For n large, this will hardly change the results and we will therefore 
ignore this effect. Now connect one of the half-edges uniformly at random to one of the remaining L n — 1 
half-edges. Repeat this procedure until all half-edges have been connected. 

Dembo and Montanari studied a slightly different version of the configuration model in [8], where 
the degrees of the vertices are deterministic instead of random. They prove that in that case, when the 
empirical degree distribution has finite mean, the graph sequence is also locally tree-like and uniformly 
sparse. Their proof can easily be adapted to show that this also holds for the above version of the 
configuration model using the strong law of large numbers. 

In [8j, it was also shown that the Erdos-Renyi random graph is uniformly sparse and is locally tree-like 
with an asymptotic degree distribution that is Poisson distributed (and note that for a Poisson distribution 
P = p). Therefore, clearly, the Erdos-Renyi random graph has a finite mean degree distribution. 

That these results hold for a wide variety of random graph models is not a surprise. It is believed that 
the behavior of networks shows a great universality. Distances in random graph models, for example, also 
show a remarkably universal behavior. See, e.g., [13] for an overview of results on distances in power-law 
random graphs. These distances mainly depend on the power-law exponent and not on other details of the 
graph. Note, however, that the results above only apply to graphs that converge locally to a homogeneous 
tree and thus, for instance, not for many inhomogeneous random graphs studied in [5] where the local 
structure is a multi-type Galton- Watson branching process instead. Certain parts of our proof easily 
extend to this case. 

In this paper we study smooth observables, we defer the investigation of the critical nature to a later 
paper. There we will study the behavior around the critical value /3 C where certain quantities (e.g. the 
susceptibility) have singularities. Of special interest is the critical behavior when r E (2, 3), where f3 c = 0, 
i.e., where the system is always in the ferromagnetic regime for any finite temperature. 

1.4 Overview of the proof and organization of the paper 

In this section, we give an overview of the proof of Theorem 11.41 and reduce it to the proofs of Propo- 
sitions 11.71 11.81 and 11.91 below. Proposition 11.71 establishes that the recursive relation that gives the field 
acting on the root of the infinite tree T{P, p, oo) is well-defined, in the sense that the recursion admits a 
unique fixed point h*. Proposition 11.81 is instrumental to control the implicit dependence of the pressure 
of the random Bethe lattice (f((3, B) on the inverse temperature /3 via the field h*. This is used in Propo- 
sition [L9] which proves that the derivative of the pressure with respect to j3, namely minus the internal 
energy, converges in the thermodynamic limit to the derivative of <p(/3,B). We also clearly indicate how 



our proof deviates from that by Dembo and Montanari in |7j. 

We will first analyze the case where B > and deal with B < later. We start by investigating the 
distributional recursion (|1.13j) : 

Proposition 1.7 (Tree recursion). Let B > and let {Kt}t>\ be i-i.d. according to some distribution p 
and assume that K\ < oo, a.s. Consider the sequence of random variables {h^'}t>o defined by h^> = B 
and, fort > 0, by (|1.13j) . Then, the distributions of hS 1 ' are stochastically monotone and M*' converges 
in distribution to the unique fixed point h* of the recursion (11.130 that is supported on [0, oo). 

We can now investigate the thermodynamic limit of the pressure. By the fundamental theorem of 
calculus, 



lim ip n (f3,B) = lim 



i; n (0,B)+J A^'^d/T 

M0,B) + J —^ n ((3>,B)df3' + J —^ n (f3',B)d/3' , (1.24) 

for any < e < /3. For all n > 1, we have that 

^(0, B) = log(2 cosh(5)) = 99(0, B), (1.25) 



lim 

n— >oo 



so this is also true for n — > 00. 

By the uniform sparsity of {G n } n >i, 







T ^V,B) 



- E 



cr,-cr 



»"J, 



fin 



{hi)&E n 



< 



\E„ 



n 



<c, 



(1.26) 



for some constant c. Thus, uniformly in n, 



d_ 
d~P 



7 M/3',B)d/3' 



<ce. 



(1.27) 



Using the boundedness of the derivative for /3' E [e, /3], we also have that 



pQ o /»/3 o 

lim/ —^ n (p',B)df3'= lim —MP',B)d(3'. 



(1.28) 



For /3 > 0, we will show that the partial derivative with respect to /3 of ifi n (P, B) converges to the partial 
derivative with respect to j3 of tp(j3,B). For this, we need that we can in fact ignore the dependence of 
h* on /3 when computing the latter derivative as we shall show first: 

Proposition 1.8 (Dependence of (p on (i3,B) via h*). Assume that the distribution P has strongly finite 
mean. Fix B\,B2 > and < i3\ , fi% < 00 . Let h\ and h\ be the fixed points of (J1.13P for (/3i , B\ ) and 
(02, B 2 ), respectively. Let iph*(/3,B) be defined as in (|1.12|) with {/ij}j>i replaced by i.i.d. copies of the 
specified h* . Then, 

(a) For B\ = B%, there exists a \\ < 00 such that 



\<p h *(lh,B 1 )-<p h *(p 1 ,B 1 )\<\ 1 \p 1 -fc 



IT-l 



(1.29) 



(b) For (3i =02, there exists a A2 < 00 such that 



\(Ph*(A^i)-Vhi(Pi,Bi)\ < \2\B1-B. 



IT-1 



(1.30) 



Note that this proposition only holds if r G (2,3). For r > 3, the exponent r — 1 can be improved 
to 2, as is shown in [7], but this is not of importance to the proof. We need part (b) of the proposition 
above later in the proof of Corollary 11.61 

Proposition 1.9 (Convergence of the internal energy). Assume that the graph sequence {G n } n >i is 
locally tree-like with asymptotic degree distribution P, where P has strongly finite mean, and is uniformly 
sparse. Let f3 > 0. Then, a.s., 

lim ^ n (p,B) = —cp(/3,B), 

n-Kx> dp Op 



(1.31) 



where ip(/3,B) is given in (|1.12p . 

By Proposition 11.91 and bounded convergence, 



f*P £) PP £) 

/ lim — MP',B)dp'= / —<p([3 f ,B)dp' = cp(P,B)-ip(e,B), (1.32) 

again by the fundamental theorem of calculus. 

Observing that < tanh(/i*) < 1, one can show that, by dominated convergence, (p(/3,B) is right- 
continuous in (3 = 0. Thus, letting e | 0, 



lim ip n ({3,B) = lim lim 

n— i-oo e^O n— >oo 



M0,B)+ J —MP\BW+ ip n (f3>,B)df3' 



dp' 



dp' 



<p(0, B) + lim {<p(p, B) - <p{e, B)) = <p(p, B), 

£\0 



(1.33) 



which completes the proof for B > 0. 

The Ising model with B < is equivalent to the case B > 0, because one can multiply all spin variables 
{ a i}i&[n} anci B with —1 without changing Boltzmann distribution f 1 1 . lj) . Furthermore, note that, 



^MP,B) 



n ^ 

iGfnl 



(7; 



/In 



< 1, 



so that B 1— >■ ip n {j3,B) is uniformly Lipschitz continuous with Lipschitz constant one. Therefore, 



lim Vn(/3,0) 



lim lim -;/;„( /3,i?) 
n-»oo BJ,0 



lim lim rbJB.B) =\imip(8,B). 



(1.34) 



(1.35) 



D 



The proof given above follows the line of argument in [7], but in order to prove Propositions 11.71 11.8 
and 11.91 we have to make substantial changes to generalize the proof to the infinite variance case. 

To prove Proposition 11.71 we adapt the proof of Dembo and Montanari by taking the actual forward 
degrees into account, instead of using Jensen's inequality to replace them by expected forward degrees, 
which are potentially infinite. This also makes a separate analysis of nodes that have zero offspring 
superfluous, which considerably simplifies the analysis. 

The proof of Proposition 11.81 (a) is somewhat more elaborate, because we have to distinguish between 
the cases where L in (I1.12p is small or large, but the techniques used remain similar. By, again, taking 
into account the actual degrees more precisely, the analysis is simplified however: we, for example, do not 



rely on the exponential decay of the correlations. Part (b) of this proposition is new and can be proved 
with similar techniques. The proof of Proposition 11.91 is proven in a similar way as in [7J. 

The remainder of this paper is organized as follows. First we shall review some preliminaries on 
Ising models in Section [2j Next, we shall study the tree recursion of (|1 . 13[) and prove Proposition 11.71 
in Section [3] and Proposition 11.81 in Section [H Finally, in Section [5j we shall prove Proposition 11.91 In 
Section [6] we shall study the thermodynamic quantities to prove Corollary 11.61 

2 Preliminaries 

The first result on ferromagnetic Ising models we will heavily rely on is the Griffiths, Kelly, Sherman 
(GKS) inequality, which gives various monotonicity properties: 

Lemma 2.1 (GKS inequality). Consider two Ising measures fi and [i! on graphs G = (V,E) and G = 
iy,E'), with inverse temperatures (3 and f3' and external fields B_ and B_' , respectively. If E C E' , /3 < j3' 
and < Bi < B[ for all i &V, then, for any U C V, 

A weaker version of this inequality was first proved by Griffiths [T2] and later generalized by Kelly and 
Sherman |15j . The second result on ferromagnetic Ising models is an inequality by Griffiths, Hurst and 
Sherman [13] which shows the concavity of the magnetization in the external (positive) magnetic fields. 

Lemma 2.2 (GHS inequality). Let (3 > and Bi > for all i € V. Denote by 

mj(B) = ix({a : a-j = +1}) - /x({<r : Oj = -1}) (2.2) 

the magnetization of vertex j when the external fields at the vertices are B_. Then, for any three vertices 

j, k,leV, 

The final preliminary observation we need is a lemma that reduces the computation of the Ising 
measure on a tree to the computation of Ising measures on subtrees: 

Lemma 2.3 (Pruning trees). For U a subtree of a finite tree T, let dU be the subset of vertices ofU that 
connect to a vertex in W = T\U. Denote by {o~ u ) the magnetization of vertex u € dU of the Ising 

model on W U {u}. Then, the marginal Ising measure on U, fijj, is the same as the Ising measure on U 
with magnetic fields 

f atanh((a u ) w J, u e dU, 

u \ B u , u£U\OU. { ' ' 

The proof of this lemma follows from a direct application of the Boltzmann distribution given in (jl.lj) , 
see [7J Lemma 4.1]. 



3 Tree recursion: proof of Proposition 11.7 



To prove Proposition 11.71 we wm first study the Ising model on a tree with I generations, T(£), with 
either + or free boundary conditions, where the Ising models on the tree T(£) with +/free boundary 



conditions are defined by the Boltzmann distributions 

/^ ,+ (^) = z e,+(o B \ ex P \ P Yl < J i (J j+ Yl Bi(Ti f 1 K=+i, foraiuearw}, (3-1) 

' {P, -> I (i,j)GT{£) ieT(£) 



and 

" ^ V, (3.2 



^ /(fT) = W^l) exp ^ 



(iJ)€T(£) ieT(£) 

respectively, where Z e ' + 'f are the proper normalization factors and dT(£) denotes the vertices in the 
£-th generation of T(£). In the next lemma we will show that the effect of these boundary conditions 
vanishes when £ — > oo. This lemma is a generalization of Lemma 4.3], where this result was proved in 
expectation for graphs with a finite-variance degree distribution. This generalization is possible by taking 
the degrees into account more precisely, instead of using Jensen's inequality to replace them by average 
degrees. This also simplifies the proof. 

We will then show that the recursion (|1.13p has a fixed point and use a coupling with the root 
magnetization in trees and Lemma 13.11 to show that this fixed point does not depend on the initial 
distribution hS°\ thus showing that (|1.13p has a unique fixed point. 

Lemma 3.1 (Vanishing effect of boundary conditions). Let m ,+ >->(B_) denote the root magnetization 
given T{£) with external field per vertex B; L > B m \ a > when the tree has +/free boundary conditions. 
Assume that the forward degrees satisfy Aj < oo a.s., for all i E T(£ — 1). Let < f3 < /3 max < oo. Then, 
there exists an M = M(/3 max , -B m j n ) < oo such that, a.s., 



m e >+(B) - m ej (B) < ^-, (3.3) 



for all £> 1. 



Remark. Lemma \3.1\ is extremely general. For example, it also applies to trees arising from multitype 
branching processes. 

Proof. The lemma clearly holds for /3 = 0, so we assume that /3 > in the remainder of the proof. 

Denote by m (B, H) the root magnetization given T(£) with free boundary conditions, when the 
external field on the vertices i G dT(£) is Bi + Hi and Bi on all other vertices i E T(£ — 1). Condition on 
the tree T{£) and assume that the tree T{£) is finite, which is true a.s., so that we can use Lemma 12.31 
Thus, for 1 < k < £, 

m k > + (B) = m k (B, oo) = m k ~ 1 (B, {/3AJ), (3.4) 

where Aj is the forward degree of vertex i E dT(k — 1). By the GKS inequality 

m* _1 QL{0Ai}) < m k ~ 1 (B,oo). (3.5) 

Since the magnetic field at all vertices in dT{k) is at least B m i n we can write, using Lemma 12.31 and 
the GKS inequality, that 

m k ' f (B) = m k (B,0) > m k - 1 (B,^{A i }), (3.6) 

where 

£ = £(/3, -Bmin) = atanh(tanh(/3) tanh(i? m ; n )). (3.7) 

This inequality holds with equality when Bi = B m \ n for all i E dT(k). Using the GKS inequality again, 
we have that 

m k - 1 (B^{A i }) > m k - x (B,<S). (3.8) 
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Note that < £(0,B m m) < /?• Since H h-» m k (B_, H{ A j}) is concave in ff because of the GHS 
inequality, we have that 

m^GB,/?!^}) -m k -\B,0) < M (m^ 1 (B, £{\}) ~ m h -\B,0)) , (3.9) 

where 

M = M(p max ,B min )= sup —^4 V <0 °- ( 3 - 10 ) 

0</3</3max UPi-DminJ 

Thus, we can rewrite m k,+ (B) using (|3.4p and bound m k, f(B) using (|3.6p and (|3.8p . to obtain 

m k - + {B) - m k >f(B) < m k - 1 (B,/3{A i }) - m k - 1 (B,0). (3.11) 

By (J3.9D . we then have that 

m fe '+( : B) - m kJ (B) < M (m k - 1 (B^{A i }) - m k - X (B,0)\ < M (m k (B,0) - m k - 1 {B,0)) , (3.12) 

where we have used (|3.6|) in the last inequality. 

By (|3.4p and (|3.5p . m ,+ (B) is non-increasing in k and, by (|3.6p and (|3.8p . m "(B) is non-decreasing 
in A:. Thus, by summing the inequality in (13.12P over k, we get that 

t e 

£(m e ' + (B)-m ej (B)'\ < J^ (m k ' + (B) - m k ' f (B_j) < M^ (m k (B, 0) - m k ~\B, 0)' 

fc=i fc=i 

= Af (m £ (S, 0) - m°CB, 0)) < M, (3.13) 



since < m l /°(B, 0) < 1. D 

We are now ready to prove Proposition 11.71 

Proof of Proposition \1.7[ Condition on the tree T(p, oo). Then h"' = atanh(m t '-' (B)) satisfies the recur- 
sive distribution (|1.13[) because of Lemma [2T3l Since, by the GKS inequality, m*'^(i3), and hence also 
h"\ are monotonically increasing in t, we have that B = h^- ' < hP' < B + Dq < oo for all t > 0, where 
Dq is the degree of the root. So, h^> converges to some limit h. Since this holds a.s. for any tree T(p, oo), 
the distribution of h also exists and one can show that this limit is a fixed point of (jl.l3|) (see [7, Proof 
of Lemma 2.3]). 

In a similar way, h^ t,+ > = atanh(m*' + (i3)) satisfies (|1.13|) when starting with h^ 0,+ ' = oo. Then, h^' + ' 
is monotonically decreasing and, for t > 1, B < hS 1 ' < B + Dq < oo, so hS t,Jr > also converges to some limit 
h. 

Let hbe a fixed point of I11.13J) . condition on this h and let h^ '*' = h. Then hS 1 '*' converges as above 
to a limit h* say, when applying f 1 1 . 1 3 1) . Note that bP' < h^'*> < h^ 0,+ K Coupling so as to have the 
same {K t }t>i while applying the recursion (11. 13H . this order is preserved by the GKS inequality, so that 
fr(*) < /»(*•*) < &(*>+) for all t > 0. By Lemma EH 



|tanh(/iW) - tanh(/i( i '+))| = \m tJ (B) - m^ + (B)\ -)• 0, for t -> oo. (3.14) 

Since the above holds a.s. for any tree T(p, oo) and any realization of h* , the distributions of h, h and h* 
are equal, and, since h is a fixed point of (|1.13p . are all equal in distribution to h. D 
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4 Dependence of cp on ({3, B) via h*: proof of Proposition 11.81 



We will now prove Proposition 11.81 by first bounding the dependence of (p on h* in Lemma 14.11 and 
subsequently bounding the dependence of h* on p and B in Lemmas 14.21 and 14.31 respectively. 

Lemma 4.1 (Dependence of (p on h*). Assume that distribution P has strongly finite mean. Fix Bi,B 2 > 
and < Pi, fa < oo. Let h* and h\ be the fixed points of (|1.13p for (Pi,B{) and (fa,B2), respectively. 
Let (ph*({3,B) be defined as in (|1.12|) with {/ij}j>i replaced by i.i.d. copies of the specified h* . Then, for 
some A < oo, 

\ip h *(Pi,Bi)-ip h *(fa,Bi)\ < Alltanh^^-tanh^)^" 1 , (4.1) 

where \\X — Y\\ MK denotes the Monge-Kantorovich-Wasserstein distance between random variables X and 
Y , i.e., \X — 1"||mk is the infimum ofE[\X — Y\] over all couplings (X, Y) of X and Y. 

Proof. Let Xi and Yi be i.i.d. copies of X = tanh(/ij) and Y = tanh^/iJj) respectively and also independent 
of L. When \\X — Y\\ MK = or \\X — y|| MK = oo, the statement in the lemma clearly holds. Thus, without 
loss of generality, we fix 7 > 1 and assume that (Xi, Yi) are i.i.d. pairs, independent of L, that are coupled 
in such a way that E|Xj — Yi\< j\\X — V||mk < °°- 
Let p = tanh(/3i) and, for £ > 2, 



F e {x 1 ,...,x e )=logle B fl(l + f3x t ) + e- B fl(l-f3x l )\-j^j ^ Iog(l 4-^^), 

I i=l i=l J l<i<j<? 



(4.2) 



and let 



Fi(xi,x 2 ) = -{log(e B (l + /3xi) + e- B (l-/3xi)) 



(4.3) 



+ log (e B (l + Px 2 ) + e' B (l - Px 2 )) - log(l + px x x 2 ) ' 

Then, with L having distribution P , 

<p h *(Pl,Bi) = F +E[F L (Xi,...,X max{2iL} )} and ^(A, B x ) = F + E[F L (Y lt . . . ,y max{2 , L} )], (4.4) 

for some constant Fq that is independent of p and B. In the remainder of the proof we will assume that 
F\ is defined as in (|4.2p . The proof, however, also works for F\ as defined in f|4.3j) . 

We will split the absolute difference between iph*(Pi,Bi) and iph*(Pi,Bi) into two parts depending 
on whether L is small or large, i.e., for some constant 9 > to be chosen later on, we split 



E[F L (Yi, ...,Y L )- F L (Xi, . . .,X L )] < E[(F L (Y U . . . ,Y L ) - F L (X h . . .,X L ))t 



[ {L<9}\ 



+ 



E[(F L (Yi,...,Y L )-F L (Xi,...,X L ))t {L>e} ] . (4.5) 



Note that 
F £ (Yi,.. 



l Y i )-F e (Xi 



,-X>)= / ^-F e (sYi + (l-s)Xi,...,sY e + (l-s)X e 
Jo ds 

= / V(y 4 -x i )-^(tyi + (i-t)Xi,...,ty / + (i-t)x < )dt 

JO i=l ax i 

= J2( Y i- x i)l -^(tYi + (l-t)Xi,...,tY e + (l-t)X e )dt. (4.6) 
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As observed in [TJ Corollary 6.3], -^ is uniformly bounded, so that 



F t {y u ...,Y e )- F t {X u ... ,X t )\ < X 1 J2 \Yi ~ Xi\, 
where Ai is allowed to change from line to line. Hence, 



(4.7) 



i=l 



E[(F L (Yi, . . . , Y L ) - F L (Xi, . . .,X L ))1 



{L>e}\ 



< E 



^\Yi- Xi\ Cl \ {L>e} 



i=l 



<Ai||X-Y|| MK E[Ll {L> } ]. 



(4.8) 



We compute, using that L > 0, 



0+1 



E[Lt {L>e} ] = ^P[Ll {L>e} >x} = ^P[Ll {L>e} > x] + Y; nLt { L>e} > x] 

x=l x=l x=6+2 

6+1 oo oo 

= J>[L>0 + 1] + ^ P[L > x] < {9 + 1) • c(0 + l)"^- 1 ) + Yj cx-^-V 
x=i x=e+2 x=e+2 



so that 



E[(F L (Y 1 , ...,Y L )- F L (X 1 , ..., X L ))t {L>e} ] 
By the fundamental theorem of calculus, we can also write 



< XAX-Y] 



-(r-2) 



F t (Yi, ...,Yi)- F l {X 1 , ...,X i ) = J2 Aii^i + Y,( Y i ~ X *)( Y J ~ X j)fij 



W 



i=i 



with 



and 



AiF t = (Xi - Xi) / -^{X x , . ..,tY { + (l- t)X i} . . .,X e )dt, 

JO OXi 



./; 



(0 



i rt pV 



&F, 



13 Jo Jo dxidxj 
Therefore, 



E[(F L (Y U ...,Yl)~ F L (X h . . . ,X L ))t {L < e} ] 



< 



E 



+ 



J2 A i F L^{L<e} 

L 

YiYi-XiWj-Xj)/^!^ 



E 



i^i 



Since 9x .q x . is also uniformly bounded ([71 Corollary 6.3]), we obtain 



E 






<\ 2 E ^2\Y i -X i \\Y j -X j \l {L < ie} 
<X 2 \\X-Y\\ 2 MK E[L 2 t {Lm ], 



(4.9) 
(4.10) 

(4.11) 

(4.12) 



[sYx + (1 - s)X 1 , ...,sYi + (l- s)Xi, ...,sY e + (l- s)X i )dsdt. (4.13) 



(4.14) 



(4.15) 
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where A2 is allowed to change from line to line. The second moment of a non-negative integer-valued 
random variable M, can be written as 



E[M 2 ] = Y^(2x - 1)P[M > 



x 



(4.16) 



x=l 



so that 



E[L 2 l {L < e} ] = J2( 2x ~ lW[Lt {L <e} >x} = Y^{2x - 1)P[L1 {L < 0} > x] 

x=l x=l 

e e 

< Y^ 2xF[L > x] < ^2 2x ■ cx _(t_1) < A 2 0~ (t ~ 3) . 



x=l 



x=l 



We split 



E 



^A 4 F L 1 



{L<9} 



4 = 1 



< 



E 



E A ^ 



i=\ 



+ 



E 



^2AiF L l {L>e} 



i=l 



(4.17) 



(4.18) 



By symmetry of the functions Ft with respect to their arguments, for i.i.d. (JQ, Yi) independent of L, 



E 



E A ^ 



i=l 



E [LAiFl] = E 



L(Y 1 -X 1 ) [ ^±( t Y 1 + (l-t)X 1 ,X 2 ,...,X L )dt 

JO OXi 



(4.19) 



Differentiating (|4.2p gives, for £ > 2, 

— F £ (xi,...,^) =^(xi,^(x 2 ,...,:^)) - ^rT[E^( Xl ' x i 
where V>( x > y) = xy/(l + /3xy) and 



j/' 



(4.20) 



i=2 



ge(x 2 , . . . , x^) = tanh 1? + V^ atanh(/3i 

i=2 



■j; > 



while differentiating (|4.3[) gives 



9 

-— F^(xi,x 2 ) = i/}(xi,gi) - i/j(xi,x 2 ). 
ax\ 



Using that IPi = Ppi-i, we have that, with K distributed as /?, 

E[LiP(X 1 ,g L (X 2 , . . .,X L ))] = PE[i>{X x ,g K+l (X 2 , . . .,X K+l ))) = PE[^{X ll X 2 ) 



(4.21) 



(4.22) 



(4.23) 



because gx+i (X 2 , . . . , Xk+i) is a fixed point of (|1.13p . so that gx+i (X 2 , . . . ,Xk+i) = X 2 and is inde- 
pendent of X\. Therefore, one can show that 



dFr 

E[L-^(x,X 2 , . . . ,X max{2>L} )] = 0, for all x G [-1, 1]. 

Since -q^- is uniformly bounded, L-^- is integrable, so that, by Fubini's theorem and (|4.24|) . 



(4.24) 



E 



E A ^ 



i=l 



E 



(Y 1 -X 1 ) [ E 
Jo 



dFr 
L-^(tY 1 + {l-t)X 1 ,X 2 ,...,X L ) 



X h Yx 



At 



0. (4.25) 
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Furthermore, by (|4.12j) and the uniform boundedness of Tj-f , 



E 



^A,F L 1 { 



L>e} 



8=1 



<E 



^l^-Xildl 



{L>0} 



j=l 



<Ai||X-F| 



-(T-2) 



(4.26) 



Therefore, we conclude that 

E[(F L (Y 1 , ...,Y L )- F L (X U . . .,X L ))l {Lm ] 



< Ai||x - y|| M K0- (r - 2) + a 2 ||x - y||L^" (r_3) - (4.27) 



Combining (|4.10p and (|4.27p and letting = ||AC — y|| M K yields the desired result. 



□ 



Lemma 4.2 (Dependence of h* on j3). Fix B > and < /3i,/?2 < /?max- -^ frg an <^ ^« > where we 
made the dependence of h* on (3 explicit, be the fixed points of ()1.13|) for (j3i,B) and (/32,B), respectively. 
Then, there exists a A < oo such that 

|| tanh(/4) - tanh(/4)|| MK < A|/3i - fo\. (4.28) 

Proof. For a given tree T(p, oo) we can, as in the proof of Proposition 11.71 couple tanh(/io) to the root 
magnetizations rrig (B) such that, for all /3 > and £ > 0, 



m 



A/ 



(5) <tanh(frg) <m^* + (B), 



(4.29) 



where we made the dependence of m " ' + on /3 explicit. Without loss of generality, we assume that 
< /3i < /?2 < Anax- Then, by the GKS inequality, 



Ml 



J,+ i 



tanh(^ 2 ) - tanh^JI < m^(B) - m^{B) = m^(B) - m £ / 2 (B) + m%{B) - m%{B). (4.30) 
By Lemma 13. 11 a.s., 



m^{B)-m'£{B)< 



M 



P f 

for some M < oo. Since mi (B) is non-decreasing in (3 by the GKS inequality, 



Mi 



Mi 



m% (B) - m% (B) < (ft - ft ) sup 

/3i</3</3 max 



d(3 ' 



(4.31) 



(4.32) 



Letting £ — > oo, it thus suffices to show that dm e, f /d(3 is, a.s., bounded uniformly in £ and < (3\ < (3 < 

Pmax- 

From [7, Lemma 4.6] we know that 







m^{^B)<Y,Vk^ 



013 



k=0 



with 







V ^= £ A,— m^B,0) 



iedT(k) 



B=B 



By Lemma 12.31 and the GHS inequality, 



A m «a,o) = A m « (as) < o *- 1(ft0) , 



(4.33) 



(4.34) 



(4.35) 
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for some field H_ , so that V k £ is non-increasing in £. We may assume that B^ > B m \ n for all i € T(^) for 
some -B m i n . Thus, also using Lemma [ 



V k/ <V k , k+1 = V A,— -m fc+1 (£,0) < V A,— -m k (B,ti{&i}) 

c — ' OH; B=B i — ' ODi 

iedT(k) - iedT(fc) 



B=B 



JLm k (B,H{A t }) 



H=€(p,B min ) 



(4.36) 



where £ = £(/?, -B m in) is defined in (|3.7p . By the GHS inequality this derivative is non-increasing in if, so 
that, by Lemma |2,31 the above is at most 



m*(B,t{Ai})-m k {B,0) <- m k+i (B, 0) - m*(B,0) 



Therefore, 



„ £-1 £-i 

—m^(j5,B) < £V M <^[m W (B,0)-ffl fc (B,0) 



9/3 



fc=0 



fc=0 



< - < oo, 



for < ft < /3 < &, 



(4.37) 



(4.38) 



D 



Lemma 4.3 (Dependence of /i* on i?). Fix /3 > and B±,B 2 > B m [ n > 0. Lei /i^ and /ijL, where we 
made the dependence of h* on B explicit, be the fixed points of (|1.13p for (/3, i?i) and (f3,B 2 ), respectively. 
Then, there exists a A < oo suca £aa£ 



tanh(n^) - tanh(/t^ 2 )|| MK < A|Bi - B 2 \. 



(4.39) 



Proof. This lemma can be proved along the same lines as Lemma [4.2i Therefore, for a given tree T(p, oo), 
we can couple tanh(Zi^) to the root magnetizations m ^' + {B) such that, for all B > and £ > 0, 

m A/ (B) < t&nh(h* B ) < m £ ' + (B). (4.40) 

Without loss of generality, we assume that < B m ^ n < B\ <B 2 . Then, by the GKS inequality, 

| tanh(/i^ 2 ) - tanh(/i^)| < m e ' + (B 2 ) - m iJ {B x ) = m i ' + (B 2 ) - m ej (B 2 ) + m ej (B 2 ) - m iJ (B x ). (4.41) 

By Lemma 13. 11 a.s., 

m e '+(B 2 )-m e ' f (B 2 )<^-, (4.42) 

for some M < oo. Since m '■'(B) is non-decreasing in B by the GKS inequality, 

OrrfJ 



m lJ {B 2 ) - m lJ {B±) < (B 2 - B t ) sup — — 

B>B min >0 OB 



(4.43) 



Letting £ — >• oo, it thus suffices to show that dm ," /dB is bounded uniformly in £ and B > B mm > 0. 
This follows from the GHS inequality: 



Qrr/J dm '* 
B>B min >o dB ~ dB 



< 



B» 



R— R . -"--"mm 



m. 



>f(B min )-m e >f(B min /2) 



<— — <oo. 

-Dmin 



(4.44) 
D 
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5 Convergence of the internal energy: proof of Proposition 11.91 

We shall start by identifying the thermodynamic limit of the intensive internal energy: 

Lemma 5.1 (From graphs to trees). Assume that the graph sequence {G n } n >i is locally tree-like with 
asymptotic degree distribution P, where P has finite mean, and is uniformly sparse. Then, a.s., 



lim ^W„03,B) = 5e 

n->oo dp 2 



0-102 



(5.1) 



where v' 2 is defined in (jl.2ip . 



Lemma 15.11 shall be proved in Section 15.11 Next, we will compute the derivative of (p(/3,B) with 
respect to /3 in the following lemma and show that it equals the one on the graph: 



Lemma 5.2 (Tree analysis). Assume that distribution P has strongly finite mean. Then, 

JL (p (P,B) = ^E[(o- 1 a 2 ) v , 

where v' 2 is defined in (|1.2ip . 

Lemma 15.21 shall be proved in Section T5.21 Lemmas 15.11 and 15.21 clearly imply Proposition 11.91 

5.1 From graphs to trees: proof of Lemma 15.11 

This lemma can be proved as in [TJ. The idea is to note that 



(5.2) 



□ 



op n *-^ 



over 



\E n \ E(i,j)eE n (o^A 



^3, 



n 



\E„ 



(5.3) 



By the local convergence and the uniform sparsity, we have that, a.s. (see (jl.lOp ). 

\E„ 



lim 



n— >oo 77, 



P/2. 



(5.4) 



The second term of the right hand side of (|5.3p can be seen as the expectation with respect to a uniformly 
chosen edge (i,j) of the correlation {<x%aj) . For a uniformly chosen edge (i,j), denote by B^^{t) all 
vertices at distance from either vertex i or j at most t, and let dBuj-\(t) = Buj\(t) \Buj\(t — 1). By the 
GKS inequality, for any t > 1, 



( aia ^B (iJ) (t) ^ ^V ^ ^ a i >%,,)(<)> 



(5.5) 



where (a^aj} 

dB {i At). 



+// 



J'B, 



(i,j) 



(t) 



is the correlation in the Ising model on Buj\(t) with +/free boundary conditions on 



Let T(p,t) be the tree formed by joining the roots, (f>i and fo, of two branching processes with 
t generations and with offspring p at each vertex, also at the roots. Then, taking n — > oo, Buj\(t) 
converges to T(p,t), because of the local convergence of the graph sequence. After all, a random edge 
can be chosen, by first picking a vertex with probability proportional to its degree, and then selecting a 
neighbor uniformly at random. Using this, one can show (see [7J Lemma 6.4]), also using the uniform 
sparsity, that, for all t > 1, a.s., 



lim E 



(i,j) 



^bL® 



E 



\°>i 



'T( P ,t) 



(5.6) 
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where the first expectation is with respect to a uniformly at random chosen edge {i,j) £ E n and the 
second expectation with respect to the tree T(p,t). By Lemma 1231 and Proposition II, 7\ 



lim E 

t— >oo 



((TfaVfo) 



T(p,t) 



E 



\<*XV2, 



thus proving the lemma. 



(5.7) 
□ 



5.2 Tree analysis: proof of Lemma 15.21 

Let Xi,i > 1, be i.i.d. copies of tanh(/i*), also independent of L. Then, with L having distribution P and 
F e defined in dO) and (Ob . 



#,5) = F +E[F L (Xi, . . . ,X max{2iL} )], 



(5.8) 



for some constant Fq that is independent of j3 and 5. 

From Proposition 11.81 it follows that we can assume that j3 is fixed in h* when differentiating (p(/3, B) 
with respect to j3. Thus, taking the derivative of (|5.8j) and using (J4.23J) . one can show that 



d_ 

dp 



<p(/3,B) = ^P + ^E[^(X 1 ,X 2 )} = |e 



P + 



X\X 2 



\ + PXiX 2 \ 



p 



E 



p + x x x 2 
1 + px x x 2 



Since, with hi, ha Li.d. copies of h*, 



E 



P + X X X 2 
1 + PXiX 2 



E 



E 



tanh(/3) + tanh(/ii) tanh(/i2) 
1 + tanh(/3) tanh(/ii) tanh(/i 2 ) 
e /3+hi+h 2 _ e -0-hi+ha _ e -0+h!-h 2 + e /3-hi-h 2 



e /3+hi+h,2 _|_ e —/3—hi+h2 _|_ g—fH+hi—h-z _|_ ^p—hx—hi 

where z/ 2 is given in (|1.21|) . we have proved the lemma. 



E 



0-10-2, 



(5.9) 



(5.10) 
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6 Thermodynamic quantities: proofs of Theorem LL5] and Corollary [176] 

To prove the statements in Theorem 11.51 we need to show that we can interchange the limit of n —> oo 
and the derivatives of the finite volume pressure. We can do this using the monotonicity properties of 
the Ising model and the following lemma: 

Lemma 6.1 (Interchanging limits and derivatives). Let {f n {x)} n >i be a sequence of functions that are 
twice differentiate in x. Assume that 

(i) lim n _ ) . 0O f n {x) = f(x) for some function x h-> f(x) which is differentiable in x; 

(ii) -fafn{x) is monotone in [x — h, x + h] for all n> 1 and some h > 0. 

Then, 



lim -r-fn{x) 

n-HX) (lX 



d 

d.r 



m- 



(6.1) 
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Proof. First, suppose that -^nf n {x') > for all x' £ [x — h,x + h], all n > 1 and some h > 0. Then, for 
/i > sufficiently small and all n > 1, 



f n (x - h) - f n (x) d f n (x + h) -fnjx) 

-h ~ dx In[X) ~ h 

and taking n->oowe get, by assumption (i), that 

f{x-h)-f{x) . f d d /( a + fe)-/( a ) 
< hminf — - / n (x) < hm sup — f n {x) < . 

— h n->oo dx n->oo CiX h 

Taking h I now proves the result. The proof for ^jf n {x) < is similar. 

We are now ready to prove Theorem 11.51 
Proof of Theorem \l.h\ We apply Lemma l6.1l with the role of f n taken by B i— >■ iJ} n (/3,B), since 



(6.2) 

(6.3) 
□ 



M n (P,B) = ^Y / (a i ) IXn = ^ n (P,B) 



(6.4) 



te[n] 



and linijj^oo ip n ((3, B) = (p(/3,B) by Theorem 11.41 and -B h-> M n ((3,B) is non-decreasing by the GKS 
inequality. Therefore, 

lim Af n (/3, B) = lim ^n(P, B) = —tptf, B), (6.5) 

which proves part (a). 

Part (b) follows immediately from Proposition 11.91 and the observation that 



Un 






(i,j)eE„ 



d(3 



(6.6) 



Part (c) is proved using Lemma [6. II by combining part (a) of this theorem and that B i— > -^M n {(5, B) 
is non-increasing by the GHS inequality. □ 

We can now prove each of the statements in Corollary 11.61 by taking the proper derivative of <p(l3, B). 

Proof of Corollary li.61 It follows from Theorem 11.51 (a) that the magnetization per vertex is given by 



M(P,B) = -^<p(P,B) 



(6.7) 



Similar to the proof of Lemma [5.21 we can ignore the dependence of h* on B when differentiating <p(/3, B) 
with respect to B by Proposition 11.81 Therefore, with f3 = tanh(/3), 



^(/5,5) = Ae 
8B yKy ' ' dB 



log ( e B JJ{1 + tanh(/3) tanh(/ Ji )} + e~ B JJ{1 - tanh(/3) tanh(/i;)} 



%=\ 



i=i 



E 



E 



.b-tiL 



b-tiL 



e ij nf=i( 1 + / 3tanh (^))-e~ ij nf=i( 1 -/ 3tanh (^)) 
e S nf=i(l + /3tanh(/ l4 )) + e" s n£=i(l " P tanh(/n)) 

Ri-rL / l+/3tanh(fc t ) \ 1/2 B y,L ( l-/3tanh(fe i; ) \ 1/2 ' 

lli=l ^l-^tanhfhi)/ lli=1 Vl+/3tanh(/ii); 

ni-rL / l+/3tanh(fe t ) \ 1/2 , „_ RrT L / l-fftanh^ 1/2 



-/3tanh(hi) 



1+/8 tanh(fti) 



(6.8) 
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where L has distribution P and {/ij}i>i's are i.i.d. copies of h* , independent of L. Using that atanh(x) 
\ l°g ( i^f ) the above simplifies to 



E 



„B T~\L atanh(/3 tanh(hi)) p~^Yl^ ,,— atanh(/3tanh(/ij)) 

pBTT^ patanh(/3tanh(h;)) \ p —B T~\L atanh(/3 tanh(hi)) 



tanh B + \^ atanh(/3tanh(/ij)) 



i=l 



(6.9) 

By Lemma [2T3l this indeed equals E (o"o) , where vl+i is given in (|1.18|) . which proves part (a). 

Part (b) immediately follows from Theorem 11.5( b) and Lemma 15.21 □ 

Acknowledgements. The work of RvdH and SD is supported in part by The Netherlands Organisation 
for Scientific Research (NWO). CG has been working at the Eindhoven University of Technology (TU/e) 
during the execution of this work; he acknowledges the support of the Department of Mathematics and 
Computer Science of the TU/e. 

References 

[1] R. Albert and A.-L. Barabasi. Statistical mechanics of complex networks. Reviews of Modern Physics, 
74:47-97, (2002). 

[2] D. Aldous and J.M. Steele. The objective method: probabilistic combinatorial optimization. In: 
Probability on Discrete Structures, pages 1-72, Springer, Berlin, (2004). 

[3] R.J. Baxter. Exactly Solved Models in Statistical Mechanics. Academic Press, London, (1982). 

[4] S. Bhamidi, R. van der Hofstad, G. Hooghiemstra. First passage percolation on random graphs with 
finite mean degrees. The Annals of Applied Probability, 20(5):1907-1965, (2010). 

[5] B. Bollobas, S. Janson and O. Riordan. The phase transition in inhomogeneous random graphs. 
Random Structures and Algorithms, 31:3-122, (2007). 

[6] S. Chatterjee and R. Durrett. Contact processes on random graphs with power law degree distribu- 
tions have critical value 0. Annals of Probability, 37(6):2332-2356, (2009). 

[7] A. Dembo and A. Montanari. Ising models on locally tree-like graphs. The Annals of Applied 
Probability, 20(2):565-592, (2010). 

[8] A. Dembo and A. Montanari. Gibbs measures and phase transitions on sparse random graphs. 
Brazilian Journal of Probability and Statistics, 24(2):137-211, (2010). 

[9] S.N. Dorogovtsev, A.V. Goltsev and J.F.F. Mendes. Ising models on networks with an arbitrary 
distribution of connections. Physical Review E, 66:016104, (2002). 

[10] S.N. Dorogovtsev, A.V. Goltsev and J.F.F. Mendes. Critical phenomena in complex networks. Re- 
views of Modern Physics, 80(4):1275-1335, (2008). 

[11] W. Feller. An Introduction to Probability Theory and Its Applications, Volume II, 2nd Edition. John 
Wiley and Sons, New York, (1971). 

[12] R.B. Griffiths. Correlations in Ising ferromagnets. II. External magnetic fields. Journal of Mathe- 
matical Physics, 8(3):484-489, (1967). 

20 



[13] R.B. Griffiths, C.A. Hurst and S. Sherman. Concavity of magnetization of an Ising ferromagnet in 
a positive external field. Journal of Mathematical Physics, 11(3):790— 795, (1970). 

[14] R. van der Hofstad and G. Hooghiemstra. Universality for distances in power-law random graphs. 
Journal of Mathematical Physics, 49:125209, (2008). 

[15] D.G. Kelly and S. Sherman. General Griffiths' inequalities on correlations in Ising ferromagnets. 
Journal of Mathematical Physics, 9(3):466-484, (1968). 

[16] M. Leone, A. Vazquez, A. Vespignani and R. Zecchina. Ferromagnetic ordering in graphs with 
arbitrary degree distribution. The European Physical Journal B, 28:191-197, (2002). 

[17] A. Montanari, E. Mossel and A. Sly. The weak limit of Ising models on locally tree-like graphs. To 
appear in: Probability Theory and Related Fields, (2010). 

[18] M.E.J. Newman. The structure and function of complex networks. SIAM Review, 45(2):167-256, 
(2003). 

[19] M. Niss. History of the Lenz-Ising model 1920-1950: from ferromagnetic to cooperative phenomena. 
Archive for History of Exact Sciences, 59(3):267-318, (2005). 

[20] M. Niss. History of the Lenz-Ising Model 1950-1965: from irrelevance to relevance. Archive for 
History of Exact Sciences, 63(3):243-287, (2009). 

[21] L. De Sanctis and F. Guerra. Mean field dilute ferromagnet: high temperature and zero temperature 
behavior. Journal of Statistical Physics, 132:759-785, (2008). 

[22] S.H. Strogatz. Exploring complex networks. Nature, 410:268-276, (2001). 



21 



